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4. Matrix integrals 

Let ^jv be the space of Hermitian matrices of size N. The inner product on f)jv is given by {A, B) = 
Tr (AB). In this section we will consider integrals of the form 

where the Lebesgue measure dA is normalized by the condition / e~^''^^^^/^d^ = 1, and S{A) = 
Tr(^^)/2 — X^^xj^mTrC^™)/™ is the action functional.^ We will be interested the behavior of the 
coefficients of the expansion of Zj^ in for large A'^. The study of this behavior will lead us to considering 
not simply Feynman graphs, but actually fat (or ribbon) graphs, which are in fact 2-dimensional surfaces. 
Thus, before we proceed further, we need to do some 2-dimensional combinatorial topology. 

4.1. Fat graphs. Recall from the proof of Feynman's theorem that given a finite collection of flowers 
and a pairing a on the set T of endpoints of their edges, we can obtain a graph Fo- by connecting (or 
gluing) the points which fall in the same pair. 

Now, given an i-flower, let us inscribe it in a closed disk D (so that the ends of the edges are on 
the boundary) and take its small tubular neighborhood in D. This produces a region with piecewise 
smooth boundary. We will equip this region with an orientation, and call it a fat i-valent flower. The 
boundary of a fat i-valent flower has the form P1Q1P2Q2 ■ ■ ■ PiQiPi, where Pi, Qi are the angle points, 
the intervals PjQj are arcs on dD, and QjPj+i are (smooth) arcs lying inside D (see Fig. 10). 




Now, given a collection of usual flowers and a pairing a as above, we can consider the corresponding 
fat flowers, and glue them (respecting the orientation) along intervals PjQj according to a. This will 
produce a compact oriented surface with boundary (the boundary is glued from intervals PjQj+i). 

We will denote this surface by T^, and call it the fattening of F with respect to a. A fattening of 
a graph will be called a fat ( or ribbon) graph. Thus, a fat graph is not just an oriented surface with 
boundary, but such a surface together with a partition of this surface into fat flowers. 

Note that the same graph F can have many different fattenings, and in particular the genus g of the 
fattening is not determined by F (see Fig. 11). 

4.2. Matrix integrals in large limit and plcuietr graphs. Let us now return to the study of the 
integral Z^. By the proof of Feynman's theorem, 

_n, fcri,(j/2-l) 
n *■ <T 

where the summation is taken over all pairings of T = r(n) that produce a connected graph and 
F„ denotes the contraction of the tensors Tr(A') using a. 

For a surface S with boundary, let z/(S) denote the number of connected components of the boundary. 

Proposition 4.1. = N^^^'h 



'Note that we divide by m and not by m!. We will see below why such normalization will be more convenient. 
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Ts 9 = 1 

Figure 11. Gluing a fat graph from fat flowers 

Proof. Let e, be the standard basis of C^, and e* the dual basis. Then the tensor Tr(A'") can be 
written as 

N 

One can visualize each monomial in this sum as a labeling of the angle points Pi, Qi, . . . , P^, Qm on 
the boundary of a fat m-valent flower by zi, 12, Z2, ^3, • ■ • , inn '-i- Now, the contraction using a of some 
set of such monomials is nonzero iff the subscript is constant along each boundary component of T^r 
(see Fig. 12). This imphes the result. □ 




Contraction nonzero iff 
i = j = p, j = k — r, 
k = p, i = m, 
that is 

i = r = k= p = j = m. 
Figure 12. Contraction defined by a fat graph. 

Let Gc(n) is the set of isomorphism classes of connected fat graphs with rii i-valent vertices. For 
r e G(:(n), let b{r) be the number of edges minus the number of vertices of the underlying usual 

graph r. 

Corollary 4.2. 



_ |Aut(r)| 



Proof. Let G'fat(n) = H '^'n; ^ (YlZ/iZ)^-' . This group acts on T, so that = Tg^, for any g e Giai 
(since cychc permutations of edges of a flower extend to its fattening). Moreover, the group acts 
transitively on the set of a giving a fixed fat graph Per, and the stabilizer of any a is Aut(ro.). This 
implies the result. □ 

Now for any compact surface E with boundary, let g(E) be the genus of S. Then for a connected fat 
graph f , 6(f) = 2g(f ) - 2 + j/(f ) (minus the Euler characteristic). Thus, defining Zjv(fi) = ZN{h/N), 
we find 
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Theorem 4.3. 

^ ^ ^ jY2-2g(r)fcb(r) 

This implies the following important result, due to t'Hooft. 
Theorem 4.4. (1) There exists a limit Woo '■= limjv^oo This limit is given by the formula 

reG.(n)[o] l^^^^^^l 

where G'c(n)[0] denotes the set o/ planar connected fat graphs, i.e. those which have genus zero. 
(2) Moreover, there exists an expansion In Z^/N'^ = X^g>o c^g-^"^^? where 

% = E(n^r) E ''''' 



and Gc(n)[g] denotes the set of connected fat graphs which have genus g. 

Remark 1. Genus zero fat graphs are said to be planar because the underlying usual graphs can 
be put on the 2-sphere (and hence on the plane) without self-intersections. 

Remark 2. t'Hooft's theorem may be interpreted in terms of the usual Feynman diagram expansion. 
Namely, it implies that for large N, the leading contribution to ln(Z]v(fi,/A^)) comes from the terms in 
the Feynman diagram expansion corresponding to planar graphs (i.e. those that admit an embedding 
into the 2-sphere). 

4.3. Integration over real symmetric matrices. One may also consider the matrix integral over 
the space Sn of real symmetric matrices of size N. Namely, one puts 



f^-N(N+l)/i f e-s{A)/n^^^ 

Jsn 



where S and dA are as above. Let us generalize Theorem 4.4 to this case. 

As before, consideration of the large iV limit leads to consideration of fat flowers and gluing of them. 
However, the exact nature of gluing is now somewhat different. Namely, in the Hermitian case we had 
(ej e^, Cfc e*) = SuSjk, which forced us to glue fat flowers preserving orientation. On the other hand, 
in the real symmetric case e* = e^, and the inner product of the functionals (E" ej on the space of 
symmetric matrices is given by (e^ ej,ek ei) = dik6ji + SuSjk- This means that besides the usual 
(orientation preserving) gluing of fat flowers, we now must allow gluing with a twist of the ribbon by 
180°. Fat graphs thus obtained will be called twisted fat graphs. That means, a twisted fat graph 
is a surface with boundary (possibly not orientable), together with a partition into fat flowers, and 
orientations on each of them (which may or may not match at the cuts, see Fig. 13). 




Figure 13. Twisted-fat graph 



Now one can show analogously to the Hermitian case that the 1 /N expansion of In Zm (where 
Zjv = ZN{2h/N)) is given by the same formula as before, but with summation over the set G^^(n) of 
twisted fat graphs: 
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Theorem 4.5. 

^ ^ ^ ^2-2g(r);^6(r) 

Here the genus g of a (possibly non-orientable) surface is defined by g = 1 — x/2, where x is the 
Euler characteristic. Thus the genus of RP"^ is 1/2, the genus of the Klein bottle is 1, and so on. 
In particular, we have the analog of t'Hooft's theorem. 

Theorem 4.6. (1) There exists a limit Woo := liniTv^oo ■ This limit is given by the formula 

where G*^ (n) [0] denotes the set of planar connected twisted fat graphs, i. e. those which have genus 
zero. 

(2) Moreover, there exists an expansion In Zn/N'^ = J2g>o^s^~^^' ''^here 



and G*^(n)[g] denotes the set of connected twisted fat graphs which have genus g. 

Exercise. Consider the matrix integral over the spaceqjv of quaternionic Hermitian matrices. Show 

that in this case the results are the same as in the real case, except that each twisted fat graph counts 
with a sign, equal to (—1)™, where m is the number of twistings (i.e. mismatches of orientation at 
cuts). In other words. In Zn for quaternionic matrices is equal lnZ2N for real matrices with N replaced 
by-iV. 

Hint: use that the unitary group [/(iV, H) is a real form of Sp{2N), and is a real form of the 
representation of A'^V, where V is the standard (vector) representation of Sp{2N). Compare to the 

case of real symmetric matrices, where the relevant representation is S^V for 0{N), and the case of 
complex Hermitian matrices, where it is y ® V* for GL{N). 

4.4. Application to a counting problem. Matrix integrals are so rich that even the simplest possible 
example reduces to a nontrivial counting problem. Namely, consider the matrix integral Z^ over 

complex Hermitian matrices in the case S{A) = Tr(A^)/2 — .sTr(v4^™)/2m, where = (i.e. we work 
over the ring C[s]/(s^)). In this case we can set h= 1. Then from Theorem 4.4 we get 



Tr(^2m)g-TV(A^)/2^^ ^ P™(7V), 

where Pm{N) is a polynomial, given by the formula Pm{N) = J2g>o ^s(^)^"^~^^~'^^ ' ^gi^) is the 
number of ways to glue a surface of genus g from a 2m-gon with labeled sides by gluing sides preserving 
the orientation. Indeed, in this case we have only one fat flower of valency 2m, which has to be glued 
with itself; so a direct application of our Feynman rules leads to counting ways to glue a surface of a 
given genus from a polygon. 

The value of this integral is given by the following non-trivial theorem. 

Theorem 4.7. (Harer-Zagier, 1986) 

(2m)! -A /m\ ^r,x{x - l)...{x-p) 



(P + I)! 



p=0 

The theorem is proved in the next subsections. 
Looking at the leading coefRcient of P^, we get. 

Corollary 4.8. The number of ways to glue a sphere of a 2m-gon is the Catalan number Cm = 

1 / 2m,\ 



m+l 



m 
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Corollary 4.8 actually has another (elementary combinatorial) proof, which is as follows. For each 

pairing a on the set of sides of the 2TO-gon. let us connect the midpoints of the sides that are paired by 
straight hues (Fig. 14). It is geometrically evident that if these lines don't intersect then the gluing will 
give a sphere. We claim that the converse is true as well. Indeed, we can assume that the statement is 
known for the 2m — 2-gon. Let cr be a gluing of the 2TO-gon that gives a sphere. If there is a connection 
between two adjacent sides, we may glue them and go from a 2m-gon to a 2m — 2-gon (Fig. 15). Thus, 
it is sufficient to consider the case when adjacent sides are never connected. Then there exist adjacent 
sides a and h whose lines (connecting them to some c, d) intersect with each other. Let us now replace 
(7 by another pairing cr', whose only difference from a is that a is connected to h and c to d (Fig. 16). 
One sees by inspection (check it!) that this does not decrease the number of boundary components of 
the resulting surface. Therefore, since g gives a sphere, so does cr'. But a' has adjacent sides connected, 
the case considered before, hence the claim. 




Figure 14. Pairing of sides of a 6-gon. 



1 




4 

Figure 15 




Figure 16 



Now it remains to count the number of ways to connect midpoints of sides with lines without 
intersec;tions. Suppose we draw one such line, such that the number of sides on the left of it is 2k and 
on the right is 21 (so that fc + / = m — 1). Then we face the problem of connecting the two sets of 2k 
and 21 sides without intersections. This shows that the number of gluings satisfies the recursion 

fc+i=m— 1 

In other words, the generating function ^ DmX™' = 1 + a; H satisfies the equation / — 1 = xf^. This 

implies that / = ^—2^^^^ which yields that D„i = Cm. We are done. 

Corollary 4.8 can be used to derive the the following fundamental result from the theory of random 
matrices, discovered by Wigner in 1955. 
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Theorem 4.9. (Wigner's semicircle law) Let f be a continuous function on K of at most polynomial 
growth at infinity. Then 

lim 4 / IV/(^/\/7V)e-'^^')/2j_ f f(^x)VI^dx. 

N^oc N 7,, ^ 27r J_2 

This theorem is called the semicircle law because it says that the graph of the density of eigenvalues of 
a large random Hermitian matrix distributed according to 'the "Gaussian unitary ensemble" (i.e. with 
density e-^^^ y^dA) is a semicircle. 

Proof. By Weierstrass uniform approximation theorem, we may assume that / is a polynomial. (Exer- 
cise: Justify this step). 

Thus, it suffices to check the result if f{x) = x^"^. In this case, by Corollary 4.8, the left hand side 
is Cm- On the other hand, an elementary computation yields ^ J^^x'^'^y/A — x^ = Cm, which implies 
the theorem. □ 

4.5. Hermite polynomials. The proof ^ of Theorem 4.7 given below uses Hermite polynomials. So 
let us recall their properties. 

Hermite's polynomials are defined by the formula 

H^{x) = {-ire^'£,e-^\ 

So the leading term of Hn{x) is (2x)". 

We collect the standard properties of Hn{x) in the following theorem. 

Theorem 4.10. (i) The generating function of Hn{x) is f{x,t) = X]„>o -f^n(a;)^ = e^^*~*^. 

(a) Hn{x) satisfy the differential equation f" — 2xf' + 2nf = 0. In other words, Hn{x)e~^^^'^ are 
eigenfunctions of the operator L = —\d'^ + ^x^ (Hamiltonian of the quantum harmonic oscillator) with 
eigenvalues n + ^. 

(Hi) Hn{x) are orthogonal: 



(iv) One has 



j= / e-''" Hm{x)H^{x)dx = 2"n!5„ 



e-^\'"^H2kix)dx = JML2^(^-™) 
, (m — ky. 



(if k > m, the answer is zero), 
(v) One has 

2V! ^2fefc!2(r-A:)! 
fc=o ^ ' 

Proof, (sketch) 

(i) Follows immediately from the fact that the operator X)(~'l)"iirS" ^^■PS a function g{x) to 
g{x-t). 

(ii) Follows from (i) and the fact that the function f{x,t) satisfies the PDE f^x — 2xfx + lift = 0. 

(iii) Follows from (i) by direct integration (one should compute Jg^f{x,t)f{x,u)e~^ dx using a shift 
of coordinate). 

(iv) By (i), one should calculate Jjj a;2"*e2^*~* e~^ dx. This integral equals 
The result is now obtained by extracting individual coefficients. 



adopted this proof from D. Jackson's notes 
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(v) By (iii), it suffices to show that 

—= / Hf.{x)H2k{x)e ^ dx- 



^ Jr ' fc!2(r-fc)! 
To prove this identity, let us integrate the product of three generating functions. By (i), we have 

^ / f{x,t)f{x,u)f{x,v)e-'''dx = e''^'^+'-+^^\ 

Extracting the coefficient of Vu^v^'^, we get the result. □ 
4.6. Proof of Theorem 4.7. We need to compute the integral 

/ Tr(A2")e-^(^')/2^A. 

To do this, we note that the integrand is invariant with respect to conjugation by unitary matrices. 
Therefore, the integral can be reduced to an integral over the eigenvalues Ai, . . . , Ajv oi A. 

More precisely, consider the spectrum map a : ()Ar — > M.^ /Sn- It is well known (due to H. Wcyl) 
that the direct image (T*fiA is given by the formula a^dA = Ce" ^. ^./^ n,<j(Ai - XjfdX, where C > 
is a normalization constant that will not be relevant to us. Thus, we have 



Pm{N) 



/E»(E.Ane-^^?/2n,<,(A.-A,)2dA 



Ve-^^?/^n,<,(Ai-A,)2dA 

To calculate the integral Jm in the numerator, we will use Hermite's polynomials. Observe that since 
Hn{x) are polynomials of degree n with highest coefficient 2", we have ni<j(Ai — Xj) = 
2-N{N-i)/2 det(iJfc(Ai)), where k runs through the set 0, 1, . . . , n - 1. Thus, we find 



Jm:=j (EAne-^^?/^n(^.-A,fdA = 
2m+NV2j^ f Aj-e-^^' HiXi - Xj^dX = 

(12) 

2m-iv(Jv-2)/2^ / Xf^e-^^' det{Hk{Xj)fdX = 

2m-N^N-2)/2^ j A2-e-^^^( {-in-W\{H.i{Xi)H,i{Xi))dX. 

(Here (— 1)*^ denotes the sign of a). 

Since Hcrmitc polynomials are orthogonal, the only terms which arc nonzero are the terms with 
a{i) = T(i) for i = 2, . . . , N. That is, the nonzero terms have cr = t. Thus, we have 



2m-iv(iv-2)/2^ f xf^e-^^H y Y\H,i{XifdX 

(■^^) N-1 -I ^oo 

2m-iV(iV-2)/2^,^^ . . . 7JV_1 ^ - / X^^'Hj {X) 



where 7, = Hi{x)'^e~^ dx are the squared norms of the Hermite polynomials. Applying this for 
m = and dividing by Jo, we find 

N-l 



A—n ^3 J-oo 



3=0 

Using Theorem 4.10 (iii) and (v), we find: 7^ = 2*i!y^, and hence 



VT^ • n i.-n 



, 2'=fc!2(j-fc)! 
j=0 fe=0 ' 
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Now, using (iv), we get 



(2m)! 



j=0 fe=0 ^ J \J J 



2™to! 



j=0 fc=0 



The sum over k can be represented as a constant term of a polynomial: 

3 



(j) =C.T.((l + ^)™(l + 2^-i)'=). 



Therefore, summation over j (using the formula for the sum of the geometric progression) yields 
(2m)! , ^,„(l + 2z-i)'^-l (2m)!,^, /to\ / N 



W^o-^C.T.((l + ^) ^^3^ )___^2f 

p=0 



We are done. 



